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Abstract. We study the consequences of superconduct- 
ing quark cores in neutron stars for the magnetic fields 
of pulsars. We find that within recent nonperturbative 
approaches to the effective quark interaction the di- 
quark condensate forms a superconductor of second kind 
whereas previously quark matter was considered as a first 
kind superconductor. In both cases the magnetic field 
which is generated in the surrounding hadronic shell of 
superfluid neutrons and superconducting protons can pen- 
etrate into the quark matter core since it is concentrated 
in proton vortex clusters where the field strength exceeds 
the critical value. Therefore the magnetic field will not 
be expelled from the superconducting quark core with the 
consequence that there is no decay of the magnetic fields of 
pulsars. Thus we conclude that the occurence of a super- 
conducting quark matter core in pulsars does not contra- 
dict the observational data which indicate that magnetic 
fields of pulsars have life times larger than lO*" years. 

Key words: dense matter - MHD - stars: interiors - 
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Recently, the possible formation of diquark con- 
densates in QCD at finite density has been rein- 
vestigated in a series of pape rs following Refs. 
( lAlford et al. 1998[ |Rapp et al. 1998|) . It has been shown 
that in chiral quark models with a nonperturba- 
tive 4-point interaction motivated from instantons 



( [Carter fc Diakonov 199*^ or nonp erturbative gluon prop- 
agators ( [Blaschke h Roberts 1998] , [Bloch et al. 1999D the 
anomalous quark pair amplitudes in the color antitriplet 
channel can be very large: of the order w 100 MeV. 
Therefore, in two-flavor QCD, one expects this diquark 
condensate to dominate the physics at densities beyond 
the deconfinement /chiral restoration transition and below 
the critical temperature 50 MeV) for the occurence 
of this "color superconductivity" (2SC) phase . In a 



three-flavor theory it has been found jAlford et al. 1999a 
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^chafer fc Wilczek 1999 ) that there can exist a color- 
flavor lo cked (CFL) phase fo r not too large strange quark 
masses ( Alford et al. 1999b ) where color superconductiv- 
ity is complete in the sense that diquark condensation pro- 
duces a gap for quarks of all three colors and flavors, which 
is of the same order of magnitude as that in the two-flavor 
case. 

The high-density phases of QCD at low temperatures 
are most relevant for the explanation of phenomena in 
rotating compact stars - pulsars. Conversely, the physi- 
cal properties of these objects (as far as they are mea- 
sured) could constrain our hypotheses about the state 
of matter at the extremes of densities. In contrast to 
the situation for the cooling behaviour of compact stars 
(Blaschke et al. 1999) where the CFL phase is dramati- 
cally different from the 2SC phase, we don't expect qual- 
itative changes of the magnetic field structure between 
these two phases. Consequently, we will restrict ourselves 
here to the discussion of the simpler two-flavor theory first. 

According to Bailin and Love (1984) the magnetic field 
of pulsars should be expelled from the superconducting 
interior of the star due to the Meissner effect and de- 
cay subsequently within 10"* years. If their arguments 
would hold in general, the observation of lifetimes for the 
magnetic field as large as lO*" years ( Makashima 1992 



Baym ct al. 1969 ) would exclude the occurence of an ex- 
tended superconducting quark matter phase in pulsars. 
For their estimate, they used a perturbative gluon prop- 
agator which yielded a very small pairing gap and they 
made the assumption of a homogeneous magnetic field. 
Since both assumptions seem not to be valid in gen- 
eral, we perform a reinvestigation of the question whether 
presently available knowledge about the lifetime of mag- 
netic fields of pulsars might contradict the occurence of a 
color superconducting phase of QCD at high densities. 

The free energy density in the superconducting quark 
matter phase with ud diquark pairing ( ~ 0+ and color 



antitriplet index p) is given by (Bailin & Love 1984) 
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where B = lotA is the magnetic induction, q the charge 
of the ud pair and the coefficients of the free energy are 
given by the following expressions 
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where t = {T — Tc) /T^. Here is the critical temperature, 
Pf the quark Fermi momentum, /i = i/pp + is the 
chemical potential (in zcroth order with respect to the 
coupling constant), dn/dE = /ipp/Tr^. 

The Ginzburg-Landau equations for relativistic super- 
conducting quarks are obtained in the usual way 



= arfp + P{d*pdp)dp + 7(-iV - qAfdp 
j = iqj{dpVd*p - d*pVdp) - 2q^-fdpd*pA . 



(3) 



In deriving the expression for the current j we have also 
used the Maxwell equation 



rots — Airj 



(4) 



The first of the Ginzburg-Landau equations @ has a so- 
lution which corresponds to the Meissner effect (Vdp — 
0, A = inside of the superconductor): 
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Fig. 1. The dependence of diquark energy gap A and 
Ginzburg-Landau parameter k on chemical potential 
/ i and density n for a NJL-type quark interaction 
( Berges fc Rajagopal 199S| ); no = 0.16 fm^^^ is the nuclear 
saturation density. The left dashed line denotes the critical 
chemical potential of the onset of quark superconductiv- 
ity (the corresponding baryon number densities are given 
on the upper scale), the right dashed line - the maximal 
values of chemical potential and density for stable stellar 
configurations. 
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For the case of weak fields [H < Hc2) one obtains the 
London equation: 



B + Aqrot rotB 







(6) 



where Aq is the penetration depth of the magnetic field 
into the superconducting quark condensate. The region of 
the change of the order parameter dp can also be deter- 
mined from (|^) via 



7C(3) 
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In the diquark condensate phase with a nonperturbative 
interaction the energy gap is A 100 MeV at /i « 400 
MeV [|, see Fig. |l|. We obtain for the coherence length 
= 0.8 X 10~^^ cm. For the penetration depth of the 
magnetic field we have 



/ — 37r/i 



3.6 X IQ-'^'^cm 



(8) 



^ within a dynamical confining quark model the di- 
quark gaps can be e ven larger than this estimate 
( Plaschke fc Roberts 199^ ). 



The thermodynamical critical field Hem that fully destroys 
the superconducting state in the case of a sup erconductor 
of the first kind is given ( Bailin fc Love 1984 ) 
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For the parameter values given above the critical field is 
Hcni ~ 8.7 X 1 0^^ G, i.e. by two or ders of magnitude larger 
than in Ref. ( |Baihn fc Love 1984| ). 



The Ginzburg-Landau parameter n which determines 
the behaviour of the s uperconductor in an external mag- 
netic field is given by ( Bailin fc Love 1984 ) 
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For values of ^ ~ pp ~ 400 MeV/c and A = 100 MeV 
we obtain k, = 34, see also Fig. 1. Therefore the supercon- 
ducting quark condensate appears as a superconductor of 
the second kind into which the external magnetic field can 
penetrate by forming quantized vortex lines in the interval 
Hci < H < Hc2 ■ The upper critical field Hc2 is determined 

by 



H'^ - 



6A2 
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The magnetic flux of the quark vortex lines <i>q amounts 
to 



27r?ic 2Trhc nhc 

= -Tr = G =6$o 

q e/S e 



(12) 



where — 2 x 10^^ G cm^ is the quantum of the proton 
magnetic flux. The lower critical field for the occurence of 
quark vortex lines is then 



67rA^ Cq 



(13) 



Here we have taken into account the spherical shape of 



the quark core in a neutron star (Sedrakian et al. 1984). 

Now we can describe the magnetic structure of a su- 
perconducting quark condensate in a pulsar and its time 
evolution. When during the cooling of the protoneutron 
star with a quark matter core the critical temperature for 
the transition to the superconducting state is reached in 
the presence of a magnetic field, then this field remains in 
the quark phase in the form of quantized vortex lines. 

At some point in the further rapid cooling of the 
star due to neutrino emission the neutrons in the 
hadronic phase ("npe" -phase) of the star become super- 
fluid. Since the basic interaction between isolated pro- 
tons resembles that of the neutrons, the protons in the 
hadronic phase become superfluid too. Since the den- 
sity of protons in "npe" -phase is only few per cent of 
the neutron density, the protons will pair in ^Sq pair- 



ing state ( Chao et al. 197^ ; Amundsen & Ostgaard 1985 



Baldo ct al. 1992). The neutrons take part in the rota- 



tion, forming a lattice of quantized vortex lines. Because 
of the strong interaction of the neutrons with the pro- 
tons a part of the super conducting protons will be en- 
trained by the ne utrons ( 3edrakian fc Shahabasian 198C ; 
Alpar et al. 1984 ) and create in the region of the neu- 
t ron vortex a magnetic field of strength H{r) given by 
( ^edrakian fc Shahabasian 198C ; Sedrakian et al. 1983 ) 



(14) 



where b = y iTh/\/3mn^ is the lattice spacing of the neu- 
tron vortex lattice, k = (m* — mp)/mp is the entrainment 
coefficient with the effective mass m* and the bare mass 
TOp of the protons; i^n is the unit vector in the direction 
of the vortex axis, r is the distance from the center of 
the vortex and H. is the angular velocity of the rotation 
of the star. This field, whose magnitude is determined by 
the rotation of the star, acts as an external field for the 
non-entrained protons and creates a cluster of proton vor- 
tices with the fluxes $o in the region around the axis of 
the neutron vortex where H{r) > H^^ . The radius of this 



region, 5^, equals ([Sedrakian et al. 1984)) 



<5„ = 6(Cp/Ap) 



(15) 



For the pulsar Vela PSR 0833-45 with = 70 rad s'^ and 
b = 10"'^ cm, we have 5^ — 10^^ cm. 

While the mean magnetic induction in the star due 
to proton vortex clusters is of the order 10^^ G, the 
mean magnetic induction within the cluster reaches 



values of 4 X 10^"* G (Sedrakian & Shahabasian 1991 



^edrakian fc Sedrakian 1995 )^ 



The magnetic field strength H{r) which occurs in the 
"npe" -phase is the strength of the external field relative 
to the superconducting quark condensate. It reaches the 
maximum value H{Q) close to the center of the neutron 
vortex, i.e. 



i7(0) = -^lnl«4.7xl0i6G 



(16) 



for Ap = 30 fm. A: = 0.7 and a coherence length = 30 
fm of the neutron. 

This external field generates quark vortex lines when 
the condition H{r) > H^^ is fulfilled. The radius of this 
region is 



5q = 5(Cq/Aq)*(^p/-^")' = 4.3 X 10"^ 



cm 



(17) 



This way, the entrainment current generates a strongly 
inhomogeneous magnetic structure in the quark conden- 
sate: the clusters of quark vortex lines with the fluxes <I>q 
and radii (5q, the axes of which are the continuation into 
the quark phase of the axes of the neutron vortex lines. 
Since 6q is by two orders of magnitude smaller than (5n, the 
mean magnetic induction in the clusters of quark vortex 
lines increases to a value of the order of 10^^ G. 

When the condition for the applicability of the London 
approximation H{Q) ^ is fulfilled, then one can apply 
the modified London equation 



B + A^rot YotB = $ai>, 



S{r 



(18) 



for the description of the magnetic structure of the quark 
condensate. The density of clusters is equal to the density 
of neutron vortex lines ny = 2 fi/Kn, where = Trh/m^ 
is the quantum of neutron circulation. 

We note that between the hadronic phase and quark 
core there is a mixed phase, in which hadrons coexist with 



a charged lattice of quark droplets (Glcndenning 1992). 
Since the number densities of neutrons and protons in the 
mixed state are lower than in the hadronic phase, these 
particles remain superfluid. So neutron vortices and clus- 
ters of proton vortex lines continue through the mixed 
phase. Therefore the magnetic field will pass through it 
and enter the quark core, see Fig. 0. In the case of small 



diquark gaps of the order of 1 MeV (Bailin fc Love 1984), 
when the diquark condensate is a superconductor of the 
first kind, the magnetic field generated in the "npe" -phase 
penetrates into the quark matter core in the form of ordi- 
nary cylindrical regions ( ^edrakian et al. 1997 ) . The radii 
of these regions will be of the order of 6q since the ther- 
modynamical critical field i?cm is of the same order as 
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Fig. 2. Magnetic field structure in the interior of a hybrid 
star with M — IAMq; b is the radius of the neutron vortex 
= 10^^ cm is the radius of the proton vortex cluster, 
and 6q = 4.3 x 10~^ cm is that of the quark vortex cluster. 
For details, see text. 



the mean magnetic field in the quark cluster. The clus- 
ters of quark vortex lines which appear due to the en- 
trainment effect in the "npe" -phase will interact with 
those which are formed by the initial magnetic field (fossil 
field). This interaction obviously implies that quark vor- 
tex lines will not be expelled from the quark core of the 
star within a time scale of r = 10^ years as suggested in 



( PaiUn fc Love 19841 ) 



We note that the evolution of the magnetic field is in- 
timately related to the rotational history of the star. In 
particular, the magnetic field of the quark core will decay 
because of the outward motion of neutron vortices when 
the star spins down. This behavior results from the fact 
that the magnetic clusters inside the quark core are the 
continuation of neutron vortices. Therefore the character- 
istic decay time of the magnetic field for the whole star 
(and also for quark core) is comparable to the pulsar's 
slowing down time, which corresponds to the life time of 
the pulsar. 

In conclusion, we find that the occurence 
superconducting quark matter core in pulsars 
not contradict the observational data which 
cate that magnetic fields of pulsars have life 
larger than 10^ 



years (Makashima 1992) 



true for small diquark gaps of the 



This 
order of 1 



of a 
does 
indi- 
times 
holds 
MeV 



(Baihn & Love 1984 


) as 


well as for 


larger ones as 


obtained recently ( 


Alford et al. 1998 




Rapp et al. 1998; 


Carter & Diakonov 1999; 


Blaschke & Roberts 199S) us- 



ing effective models for the nonperturbative quark-quark 
interaction. 



Acknowledgements. K.M.S. and D.M.S. acknowledge the hos- 
pitality of the Department of Physics at the University of 
Rostock where this research has been started. We thank H. 
Grigorian, K. Rajagopal, G. Ropke, A.D. Sedrakian and D.N. 
Voskresensky for their discussions and comments. 



, Wilczek, P., 1998, Phys. Lett. B 
Wilczek, P., 1999a, Nucl. Phys. B 



References 

Alford, M., Rajagopal, K. 

422 247 
Alford, M., Rajagopal, K., 

357 443 

Alford, M., Berges, J., Rajagopal, K., 1999b, hep-ph/9903502 
Alpar, M.A., Langer, S.A., Sauls, J.A., 1984, ApJ 282 533 
Amundsen, L., Ostgaard, E., 1985, Nucl. Phys. A 437 487 
Baihn, D., Love, A., 1984, Phys. Rep. 107 325 
Baldo, M., Cugnon, J., Lejeune, A., Lombardo, U., 1992, Nucl. 

Phys. A 536 349 
Baym, G., Pethick, C.J., Pines, D., 1969, Nature 224 673 
Berges, J., Rajagopal, K., 1999, Nucl. Phys. B 538 215 
Blaschke, D., Klahn, T., Voskresensky, D.N., 1999, astro- 

ph/9908334 

Blaschke, D., Roberts, C.D., 1998, Nucl. Phys. A 642 197 
Bloch, J.C.R., Roberts, C.D., Schmidt, S.M., 1999, Phys. Rev. 

C (in press), nucl-th/9907086 
Carter, G.W., Diakonov, D., 1999, Phys. Rev. D 60 016004 
Chao, N.-C., Clark, J., Yang, C.-H., 1972, Nucl. Phys. A 179 

320 

Glendenning, N.K., 1992, Phys. Rev. D 46 1274 
Makashima, K.: 1992, Magnetic Fields of Binary X-ray Pul- 
sars. In: The structure and evolution of neutron stars, Pines 
D., Tamagaki R., Tsuruta S. (eds.), Addison- Wesley, New 
York, p. 86 

Rapp, R., Schafer, T., Shuryak, E.V., Velkovsky, M. 

Phys. Rev. Lett. 81 53 
Schafer, T., Wilczek, F., 1999, Phys. 
Sedrakian, D.M., Shahabasian, K.M., 
Shahabasian, K.M. 



1998, 



Sedrakian, D.M., 
Afz 19 175 

Sedrakian, D.M., 
Afz 21 547 

Sedrakian, D.M., 
555 

Sedrakian, A.D., 
Sedrakian, A.D., 
RAS 290 203 



Rev. Lett. 82 3956 
1980, Afz 16 417 
Movsissian, A.G., 



1983, 



Shahabasian, K.M., Movsissian, A.G., 1984, 



Shahabasian, K.M., 1991, Sov. Phys. Usp. 34 



Sedrakian, D.M., 
Sedrakian, D.M. 



1995, Ap.J. 447 305 
Zharkov, G.F., 1997, 



MN- 



